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We present a theoretical study of the density profile of a trapped strongly interacting Fermi gas 
with unbalanced spin populations. Making the assumption of the existence of a first order phase 
transition between an unpolarized superfluid phase and a fully polarized normal phase, we show 
good agreement with a recent experiment presented in [2^ . 
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In the Bardeen-Cooper-Schrieffer (BCS) mechanism, 
the onset of superfluidity is associated to pairing between 
two fermionic species with matched Fermi levels. This 
scheme is relevant to systems like metal superconductors, 
superfluid 3 He or ultra cold gases, as observed in recent 
ground breaking experiments [jj, |2, y, Ij, 111 li| • However, 
other physical systems, like magnetized superconductors 
or neutron stars, require the understanding of fermionic 
superfluidity in the presence of mismatched chemical po- 
tentials. The nature of superfluidity in these systems 
have been the subject of longstanding debates [ilia 111; 
which have been renewed by the opportunity to address 
this topic experimentally using gaseous samples. Various 
mechanisms were proposed to describe the ground state 
of an ensemble of fermions with mismatched Fermi lev- 
els: deformed Fermi surface [13, FFLO (Fulde, Ferrell, 
Larkin, Ovshinnikov) states, where Cooper pairs acquire 
finite momentum, and their generalization to trapped 
systems 0, 0, 0, interior gap superfluidity or 
phase separation between a normal and a superfluid state 
through a first order phase transition 0, 0, 0, llij . 
When the strength of the interactions is varied, a com- 
plicated phase diag ram mixing several of these scenarios 
is expected [3111 El. 

Extending the seminal observation of fermionic super- 
fluidity in ultra-cold atom systems, two recent experi- 
ments EE have started probing the regime of mis- 
matched Fermi levels by cooling samples containing dif- 
ferent atom number in each spin state. In this paper 
we will focus on the results presented in [2^ where the 
authors studied the density profile of a gas of fermionic 
lithium when varying, in the regime of strong interac- 
tions, the population imbalance between the two trapped 
spin states. One the most striking results is displayed in 
Fig. ^ It shows that the radius of the minority com- 
ponent is strongly reduced with respect to the one of 
non-interacting gas with the same atom number. In the 
present paper, we propose an interpretation of this result 
on the basis of the existence of first order phase separa- 
tion between the normal and superfluid components and 
the use of universality in the strong interaction regime. 
We show that these two ingredients are sufficient to pro- 



vide a good quantitative agreement with experimental 
data. It is therefore complementary to recent studies of 
[3EEEEE3 j wn0 address the same topics using a more 
complicated formalism. 

Our analysis is based on the assumption of the exis- 
tence of a zero temperature first order phase transition 
between a fully polarized normal phase containing a sin- 
gle spin species and an "unpolarized" superfluid state 
composed of a balanced mixtures of the two species. The 
existence of this phase transition was suggested by pre- 
vious theoretical studies 0, fltil Il7j and is based on the 
following qualitative argument: In the grand canonical 
ensemble, the chemical potentials /x; of the two species 
are fixed. In the ground state, the system minimizes the 
grand-potential H = H — HiNi, where H is the hamil- 
tonian of the system and Ni is the population of species 
i. Let us now consider a situation where the chemical 
potentials are mismatched, with Sfj, = fJ-i — > 0. 
Promoting a particle from state 2 to state 1 decreases 
the grand potential by Sfi but implies the breaking of 
a pair, hence increases the energy by the superfluid gap 
A. When Sfi < A we therefore expect the system to be 
unpolarized, while above this threshold a fully polarized 
system is obtained. 

At zero temperature, the actual position of the phase 
transition can be found analytically in the case of a ho- 
mogeneous unitary gas by matching chemical potentials 
and pressures in the two phases |l7| . Indeed, we can 
write in the unpolarized phase ni = ni = /(/ii,/i2). 
Using the thermodynamical identity <9 Ml n2 = d^rii, we 
deduce that / can actually be expressed as a function 
of fi = (fi± + /i2)/2 only. To obtain the exact expres- 
sion for /, we then consider the case of matched Fermi 
surfaces, fi± — 112 = M- In this latter case, universal- 
ity at the unitary limit allows one to write [i = £Ep, 
where Ep — h 2 (6Tr 2 n) 2 ^ 3 /2m is the Fermi energy of a 
non-interacting gas of density n = n 12 and £ is a univer- 
sal parameter whose determination has attracted a lot of 
attention from both theoretical [rR EE EE EE| and ex- 
perimental groups 0, EE ESI EE EE| and is now thought 
to be £ ~ 0.45. Using this expression for the chemical 
potential, we then deduce that 
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Gibbs-Duhem identity dP$ = riidfii + n-id^i finally 
yields for the pressure in the superfluid phase 

In the normal phase, we assume only the majority com- 
ponent is present. We have then an ideal gas constituted 
of particles of type 1 with chemical potential /ii, thus 
giving 
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Equating P^ and Ps |3J], we see that the two phases 
coexist only if fix and satisfy the condition ^2/^1 = Vc 
with 



3/5 



-0.061, 



(4) 



a relation already found in 0,^3- ^ n a trap, the chem- 
ical potential depends on position. 

To compare with experiments, we now consider the 
case of a cloud of atoms trapped in a harmonic potential 
V{r) — m^i^iXi Without loss of generality, we 
will restrict our analysis to the case of a isotropic trap 
with frequency Q = (u x DyU z ) 1 / 3 . We can indeed always 
recover the more general anisotropic case by making the 
scaling transform Xi — > UiXi/ui. 

To calculate the density profile of the cloud, we make 
use of the local density approximation, where we assume 
that the chemical potential of species i depends on posi- 
tion as jLtj(r) — fJ-i — V(r). 

If we assume component 1 is the most populated, 
the inner superfluid region is defined by the condition 
^2{r) / m{r) < rjc and is bounded by the radius R2 de- 
fined by 
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Atoms of the minority species are located in the paired 
superfluid phase only. We thus have 
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FIG. 1: Comparison between experimental data of |23| and 
the model presented here. P = (Ni - Na)/(Ni + N 2 ) is the 
imbalance between the two spin populations. The radius R of 
each component is normalized to the radius -Rtf of the ideal 
Fermi gas containing the same atom number. Squares and 
diamonds correspond to experimental measurements of the 
radii of the minority and majority components respectively. 
The full lines are the predictions of the first order transition 
model. 
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where R = + y^/mui and 



(7) 

Excess atoms of the majority species are located be- 
tween r = R2 and r = R\ such that mu) 2 R%/2 = n\. 
The number of excess atoms is therefore N\ — N2 — 
Ir2 ni ( r ) ^ 3t *' hence 



N i~ N * = £;(ji) (9(i)-g(R 2 /Ri))- (8) 

Dividing by JHJ by © yields the implicit equation for 
r/o = lA/fJ-i as a function of N1/N2 

Ni = 1 , ,3/2 8 g(l) - gWgi) (q] 
N 2 ^ (i + vo) 3 g(R 2 /R) ' [) 

Equation @ is solved numerically and the value ob- 
tained for 770 is then used to calculate the radii R\ and 
i?2- The predicted evolution of the Ri versus the pop- 
ulation imbalance P = (iVj — N2)/(N± + N2) is shown 
in Fig. n To follow Ref. [23j, we have normalized each 
Ri to the Thomas-Fermi radius Rtf associated to an 
ideal gas containing N atoms. The agreement with the 
experimental data is quite good as soon as P > 0.1, a 
remarkable result, since the model presented here con- 
tains no adjustable parameter, as soon as the value of £ 
is known. 
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FIG. 2: Integrated column density of the majority (full line) 
and minority components (dashed line) for a population im- 
balance P = 0.57. Inset: excess column density ni — n^. By 
contrast with the experimental results, the boundary between 
the two phases is here very sharp. 



For weak population unbalance, experimental varia- 
tions of the radii is flatter than predicted by the first 
order transition model. As proposed in [23} . this sug- 
gests that the phase separation does not happen exactly 
at P = 0, but above some threshold P c ~ 0.1. This 
point is strengthened when one compare the theoreti- 
cal and experimental density profiles. As in |23j |. we 
have represented in Fig. |5]the integrated column den- 
sity rii(x) = J dym(x,y,0) for P = 0.57 (as in Fig. 
2.D of |23j). We immediately see in this figure that 
the transition between the two phases is very sharp, 
by contrast to what is observed experimentally. Finite 
temperature might explain this discrepancy. Indeed, for 
low population imbalance, the superfluid phase extends 
nearly throughout all the cloud, and in particular in re- 
gions where the density, hence the Fermi energy, is very 
low. In the unitary regime, the critical temperature for 
the superfluid-normal transition is given by the scaling 
k^T c — aEp, with a ~ 0.2 j^. It may then happen 
that near the cloud edge, the temperature T of the sam- 
ple becomes locally larger than T c . If this condition is 
satisfied at a radius R c smaller than the demixing radius 
B.2, then the first order transition to the fully unpaired 
state will not happen. If we introduce the Fermi tem- 
perature Tp — huj^Ni) 1 / 6 /ks associated to a number 
of atoms N\, the superfluid/normal and paired/unpaired 
transitions will happen at the same position if T/Tp sat- 
isfies the condition 



FIG. 3: Critical population imbalance under which finite tem- 
perature effects overule demixing. The observed experimental 
threshold P ~ 0.1 corresponds to a temperature T/Tf < 0.1 
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In fig. |3| we have plotted as a function of tempera- 
ture the evolution of the critical imbalance under which 
no demixing is expected. Experimentally, demixing only 
happens in the conditions of [2^ for P > 0.1, which, ac- 
cording to eq. corresponds to T/Tp < 0.1, a value 
compatible with experimental data. Despite a semi- 
qualitative agreement, this finite temperature argument 
needs to be clarified by a more careful analysis, following 
for instance the work presented in [5^ . Others scenarios 
can also be envisioned to explain the smooth crossover 
between the two phases, such as the existence of an in- 
termediate phase - e.g. a gapless or FFLO phases as 
proposed in [2(| -, or a breakdown of the local density 
approximation due to the fast variation of the density 
profile. Nevertheless, the good agreement between theory 
and experiment for the data presented in fig. ^ suggests 
that the crossover region between paired and unpaired 
phases should remain relatively narrow. 

From the analysis presented here, it appears that the 
observations of [23J are consistent with a scenario of a 
transition between a fully paired and a fully polarized 
phases separated by a narrow crossover region proba- 
bly related to finite temperature effects breaking the 
superfluid phase before demixing may occur. Interest- 
ingly, if this thermal scenario was confirmed the measure- 
ment of the critical population unbalance under which no 
phase separation happens could provide a useful tool for 
fermion thermometry. Another issue not addressed here 
is related to the superfluid character of the system when 
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the population imbalanced is varied. This is especially 
important if one wishes to understand the results pre- 
sented in where it is shown that a mismatch of the 
Fermi surfaces by about 50% leads to a breakdown of 
superfluidity. 
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